A T8.5 Brown Dwarf Member of the Xi Ursae Majoris System 
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ABSTRACT 

The Wide-field Infrared Survey Explorer has revealed a T8.5 brown dwarf (WISE 
J111838. 70+312537. 9) that exhibits common proper motion with a solar-neighborhood (8 pc) 
quadruple star system - Xi Ursae Majoris. The angular separation is 8.5', and the projected 
physical separation is ss 4000 AU. The sub-solar metallicity and low chromospheric activity of £ 
UMa A argue that the system has an age of at least 2 Gyr. The infrared luminosity and color of 
the brown dwarf suggests the mass of this companion ranges between 14 and 38 Mj for system 
ages of 2 and 8 Gyr respectively. 
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The effective temperature and thus spectrum of 
a brown dwarf evolves with time as it cools as a 



degenerate object (|Kumarlll962l ). For an isolated 
brown dwarf, determination of mass and age are 
intertwined, such that a broad locus of mass and 
age will be consiste nt with a single mea sured ef- 
fective temperature ([Burrows et al.ll2003l) . In rare 
instances brown dwarfs may reside in close binary 
systems, resolving this am biguity with a direct dy- 
namical mass estimation ( Konopackv et al.ll2010t 
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Cardoso et al.ll2009l ). In the absence of a dynami- 
cally measured mass, spectral modeling can signif- 
icantly constrain a brown dwarf's mass if there ex- 
ists sufficient restriction on the object's age. Con- 
straints on age and metallicity are available if the 
brown dwarf is a member of a multiple star sys- 
tem. In this case, the properties of the primary, 
particularly chromospheric activity and kinemat- 
ics, provide an indication of age. 

The Wide -field Infrared Surv ey Explorer mis- 
sion (WISE, IWright et ai1l2010l ) has been a pro- 
ductive engine for the discovery of the coolest 
brown dwarfs. The WISE Wl(3.4 /im) and W2 
(4.6 fim) filters are optimally tuned to select 
the coolest candidates, specifically those with 
spectra significantly shaped by m ethane absorp- 



tion at low effective temperature (M ainzer et al 



2005). To date spectroscopic follow-up of WISE- 
selected sources has rev ealed more than 100 ul- 
tracool brown dwarfs (Kirkp atrick et al.l 12011 ) 
including; several exceptionally cool Y-dwarfs 



( Cushing et~al1l201lllkirkpatrick et al.l2012tlTinnev et al 



20121 ) demonstrating that WISE colors provide for 
reliable photometric selection of ultra-cool brown 
dwarf candidates. 

WISE J111838. 70+312537.9, hereafter WISE 
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Table 1: Photometric Observations of WISE 
1118+31. 



Filter 



Vega Magnitude Instrument 



Fig. 1. — WISE All-sky Image Atlas cutouts show- 
ing a 10' FoV centered halfway between WISE 
1118+31 (lower left) and £ UMa (upper right). 
The left panel shows Wl (3.4 /xm), while the right 
panel shows W2 (4.6 /im). The lines point to 
WISE 1118+31. The faint source near the lines 
in the left panel is the nearby 2MASS star, #1 on 
Figure [2 The odd shape above WISE 1118+31 
in the W2 image is a ghost image of £ UMa. 
The right-most column shows 1' postage stamps 
of WISE images at 3.4, 4.6 & 12 fim. 

1118+31, easily meets t he WISE brown dwar f 
color selection criterion (jKirkpatrick et al.1 120111 ) 
with a W1-W2 color of 2.85 compared to a selec- 
tion threshold of 2.0. Few confusing objects meet 
these selection restrictions. In addition this source 
lies 8.5' from one of the nearer stars to the Sun, 
£ UMa, prompting an investigation into a possi- 
ble system membership. This paper reports the 
spectral characterization of WISE 1118+31 and 
the analysis of a series of astrometric observa- 
tions spanning 26 months aimed at determining 
whether this source exhibits common proper mo- 
tion with £ UMa. These observations demonstrate 
that WISE 1118+31 is a newly found member of 
this already remarkable multiple star system. The 
characteristics of the primary system provide an 
indication of the metallicity and age of the newly 
discovered ultra-cool brown dwarf constraining the 
mass of this object. 



2. Observations 
2.1. WISE 

WISE imaged the region of the sky containing 
WISE 1118+31 on twenty occasions between 21 
May 2010 04:20 UT and 25 May 2010 09:39 UT. 
The WISE All-Sky Catalog reports this source as 
well detected in Wl and W2, with marginal sig- 
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Note. — All magnitudes are Vega magnitudes and use 
the 2MASS JHK S filters except a s noted . The Y-band 
calibration uses the iHamuv et al.l J2006f l transformation 
of Y - K s vs. J - K s . 



nal at the position of the Wl/2 source in W3 (12 
/im) and only an upper limit in W4 (22 /im). Of 
the twenty apparitions, eighteen were sufficiently 
separated from a detector edge to permit source 
extraction. In all eighteen cases the source was 
detected in W2, producing a combined SNR=34 
detection with W2=13.31. Because exceptionally 
cool brown dwarfs are considerably fainter in Wl, 
which was optimized to produce a substantial flux 
difference between Wl and W2, WISE 1118+31 
is detected in Wl in only 12 of the 18 opportu- 
nities with a combined SNR=15 and Wl=16.16, 
yielding a color of Wl-W2=2.85. Figure [T] shows 
portions of the WISE image atlas covering both 
WISE 1118+31 and £ UMa. 

2.2. Follow-up Imaging 

Multiple epochs of near-infrared imaging pro- 
vide the astrometric data for WISE 1118+31 
needed to confirm common proper motion with 
£ UMa. Photometric information from these im- 
ages in summarized in Table [1] while the astro- 
metric data are listed in Table [2J No correc- 
tions for non-linearity were applied to the ground- 
based photometry because the sky is brighter 




Fig. 2.— Y-band image of the WISE 1118+31 
field obtained at Fan Mountain Observatory 15 
March 2012 UT. North is up and East is to the 
left. Perpendicular lines mark the position of 
WISE 1118+31. Star 1 nearby is the J=15.94 
star 2MASS J11183876+3125441, while star 2 is 
2MASS J11185086+3126520 with J=10.32. The 
faintest detected objects have Y w 20. The faint 
diffraction spike entering from the upper right 
(northwest) is from £ UMa 8.5' away. 



(jSanchez et al.l l2008t ) than both WISE 1118+31 
and the faint 2MASS objects used as calibrators. 
The abs olute calibration u ncertainties of 3% for 
Spitzer (JReach et al.l l2005h and 2.4, 2.8, 4.5 & 
5.7% for W1..4 (WISE Explanatory Supplement 
§IV.4.h.v) relative to Spitzer have not been in- 
cluded in the errors quoted in Table [TJ 

2.2.1. Fan Mountain Observatory /FanCam 

Y&J photometry and astrometry of WISE 
1118+31 were obtained at various epochs between 
28 Nov 2010 and 15 Mar 2012 with FanCam, 
a HAWAII- 1 based near- infrared imager operat- 
ing at the Univers ity of Virginia's Fan Mo untain 
31-inch telescope (JKanneganti et al.l 120091 ) . The 
source position was dithered by approximately 
10" between 30 s exposures comprising total ex- 
posure times ranging from 60 to 80 minutes. The 
FanCam field of view (FOV) is 8.7' (0.51"/pixel). 
The central 7' x 7' (Figure [2]) of the combined, 
dithered exposures was fully covered, providing 



several 2MASS stars for photometric and astro- 
metric reference. A median sky frame was sub- 
tracted from each individual exposure prior to flat 
fielding with the median background level sub- 
sequently restored to the image. Y&iJ aperture 
photometry was computed using an aperture with 
a radius of 3 pixels. The zero points for the J- 
band were computed using stars in the FOV with 
measured 2MASS magnitudes since the J-band 
filter in FanCam are based on the 2MASS system. 
In order to derive the Y-band zero point, we first 
computed the Y-band magnitudes of stars in the 
FOV using their 2MASS J a nd K„ magnitudes an d 
the transformation given by lHamuv et al.l ( 20061 ). 
The final uncertainty in the magnitudes include 
the photon noise from the sky and source, the read 
noise, and the uncertainty in the zero point due to 
the computed Y-band magnitudes of the calibra- 
tors. The resulting magnitudes and uncertainties 
are given in Table [TJ 

2.2.2. Mount Bigelow/2MASS 

The former 2MASS camera on the 1.54m 
Kuiper Telescope on Mt. Bigelow, Arizona, has 
three 256x 256-pixel NICMOS3 arrays simultane- 
ously observing in 2MASS J, H, and K s filters 
(Milligan et al. 1996). The plate scale for all 
three arrays is 1.65"/pixel, resulting in a 7' field 
of view. Exposures of 10s duration, 216 in all, 
of WISE 1118+31 were obtained on 20 May 2011 
using 6 repeats of a 3 x 3 box dither pattern, with 
four consecutive images taken at each of the nine 
dither positions. The data were reduced using 
custom IDL routines implementing standard near- 
infrared flat fielding, background removal, and co- 
addition techniques. Flat fields in each band were 
constructed using on-source frames. 2MASS stars 
provided photometric reference in all three bands, 
leading to the magnitudes reported in Table [TJ 

2.2.3. WIYN/WHIRC 

JHK S broad-band imaging of WISE 1118+31 
was obtained on UT 31 Dec 2011 with the WIYN 



High- Resolution Infrared Camera (WHIRC, Meixner et al 
l2010h and the WIYN 3.5-m Observatory. The data 
quality is excellent: both seeing (~ 0.5" with 0.1" 
pixel scale) and photometric stability conditions 
were optimal. For each band, individual frames 
had exposure times of 120, 120 and 40 seconds, for 



J, H and K S1 respectively. WHIRC uses MKO fil- 
ters for JHK S . Using an efficient on-array dither 
pattern, a total of 7, 9 and 43 frames for J, H and 
K s , respectively, were obtained, thus providing a 
total exposure time of 840, 1080 and 1720 seconds 
for the J, H and K s mosaics respectively. 

Individual frames were corrected for pupil- 
ghosts, dark and median sky flat subtracted, nor- 
malized by dome flat, and distortion corrected us- 
ing information from the WHIRC user information 
guide. Astrometric and photometric solutions us- 
ing 2MASS standards were then found. The fully 
reduced frames were combined into a deep mosaic, 
with outlier (bad pixel) rejection applied using 
temporal statistics. A final astrometric and flux 
calibration was then applied to the deep mosaic. 
The photometric uncertainty (comparing with the 
2MASS PSC) was typically better than 5% for 
each mosaic produced. The achieved spatial res- 
olution was ~ 0.5 — 0.7" for the final mosaics, 
and the astrometric uncertainty was ~ 0.05". The 
target source, WISE 1118+31, was detected in 
all three bands. Using a 1.1" radius circular aper- 
ture, the background-subtracted integrated (Vega) 
magnitudes found are reported in Table [TJ Note 
that WHIRC uses MKO filters, and the expected 
difference J2MA SS — Jmko for a T8. 5 brow n dwarf 
is « 0.36 mag ((Stephens fc Leggetd 120041 ) . If we 
correct the J2MASS values by this amount we get 
three values for Jmko: 17.79, 17.86 & 18.01. The 
mid-range of these values in Jmko = 17.9 which 
we adopt, giving Jmko ~ W2 — 4.6. 

2.24. Spitzer 



2.2.5. Keck Adaptive Optics Imaging 

High resolution imaging observations of WISE 
1118+31 were obtained using the Keck II LGS-AO 
syste m ( Wizinowich et al.l l2006t Ivan Dam et al 



The Infrared Array Camera (IRAC) (JFazio et al. 



2004) onboard the Spitzer Space Telescope em- 



ploys 256 x 256-pixel detector arrays to image a 
field of view of 5f2x5.'2 (1"2 pixel" 1 ). IRAC was 
used during the warm Spitzer mission to obtain 
deeper photometry in its 3.6 and 4.5 /im channels 
(hereafter, chl and ch2, respectively) than WISE 
was able to take in its Wl and W2 bands. These 
observations were made as part of Cycle 7 and 
Cycle 8 programs 70062 and 80109 (Kirkpatrick, 
PI). Our standard data acquisition and reduction 
methodology for IRAC o bservations is outlined in 
iKirkpatrick et al.l (|201lh . 



20061 ) with NIRC2 on 14 April 2012 (UT). WISE 
1118+3125 is not bright enough to serve as the 
tip-tilt reference star for the LGS-AO system, so 
we used the nearby 2MASS star (# 1 on Fig- 
ure [2]), which has i?=16.7 and is located about 
7" from the target. The seeing was generally 
very good throughout the night (0"3-0"5, how- 
ever high clouds were present during the observa- 
tions of WISE 1118+31. We used the MKO H 
filter and narrow plate scale (0.009942"/pixel for 
a single-frame field-of-view of 10" x 10") for the 
observations. The data were obtained by using 
a 3-point dither pattern that avoided the noisy, 
lower left quadrant of the array. Each image had 
an integration time of 120 s and the dither pattern 
was repeated five times to give a total exposure 
time of 1800 s. 

The images were reduced in a standard fashion 
using custom IDL scripts. These steps included 
dark frame subtraction, flat fielding (using a dome 
flat), and sky subtraction from a sky frame created 
from the dithered science frames. The individual 
frames were then shifted to move WISE 1118+31 
to the center of the array and the stack was me- 
dian averaged to create the final mosaic seen in 
Figure El The FWHM in the final mosaic is 54 
mas and shows no irregularities or evidence for a 
close companion. 

2.3. Spectroscopy 

2.3.1. LBT-LUCIFER 

We obtained an H- and if-band spectrum of 
WISE 1118+31 on 2012 Dec 12 (UT) using the 
Large Binocular Telescope (LBT) Near-Infrared 
Spectroscopic Utility with Camera and Integral 
Field Unit for Extr agalactic Research (LUCIFER, 
Mandel et al.ll2008l ) . A series of twelve 300 s expo- 
sures was obtained at different positions along the 
4' slit to facilitate sky subtraction. The A0 V star 
HD 97034 was also observed for telluric correction 
and flux calibration purposes. A series of halogen 
lamp exposures were also obtained for flat fielding 
purposes. 

The data were reduced using custom Interac- 



tive Data Language (IDL) softwar e based on the 



Spex tool data reduction package (jCushing et al 



2004). Pairs of images taken at two different po- 
sitions along the slit were first subtracted and 
fiat-ficldcd. The spectra were then extracted and 
wavelength calibrated using sky emission lines of 
OH and CH 4 . The twelve spectra are combined 
and then corrected for telluric absorption and 
flux calibrated usi ng the technique described by 
Vacca et al.l ( 2003 ). The final spectrum is shown 
in Figure HJ 

2.3.2. Hale - TripleSpec 

A 1-2.5 (im spectrum of WISE 1118+31 was 
obtai ned with the Trip le Spectrograph (Triple- 
spec, HerteZeFaI][200^ at the 5.08m Hale Tele- 
scope at Palomar Observatory. The 1—2.5 /xm 
range is covered over four cross-dispersed orders 
which are imaged simultaneously on the 1024 x 
2048 HAWAII-2 array. The l"-wide slit provides 
a resolving power of R ss 2700. A series of eight, 
300 s exposures were obtained at two different po- 
sitions along the 30"-long slit to facilitate sky sub- 
traction. The A0 V star HD 99966 was observed 
for telluric correction and flux calibration purposes 
and dome flats were obtained at the start of the 
night. 

The data were red uced using a modifi ed ver- 
sion of the Spextool (JCushing et al.1 120041 ) pack- 
age; a detailed description of the reduct ion steps 
can be found in iKirkpatrick et al.l ( 20111 ). Briefly, 
a two-dimensional wavelength solution is derived 
using sky emission features of OH and CH4. Spec- 
tra are then extracted from pair-subtracted, flat- 
fielded images. The resulting spectra are com- 
bined and corrected for telluric absorption and 
flux calibrated on an order-by-order basis. Finally, 
the spectra from each order are stitched together 
to form a complete 1—2.5 /xm spectrum. The 
spectrum was then flux calibrated as described in 
Ravner et al.l ( 20091 ) using the photometry in Ta- 
ble 1 . The final spectrum is shown in Figure [4] 

3. Discussion 

3.1. Properties of the Central Star System 

£ Ursae MajorisL 1 ! is a complex stellar sys- 
tem with at least 4, and possibly 5 components 



* 






Fig. 3.— ff-band LGS-AO image of WISE 
1118+3125 taken with the NIRC2 camera on Keck 
II. The image is «0'/6 on a side with North up and 
East to the left. There is no evidence for an equal 
brightness companion beyond a separation of 50 



Palomar/TripleSpec 

LBT/Lucifer 1 




ialso Alula Australis, HR 4374/5, Gl 423, HD 98230/1 



1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 
Wavelength (|xm) 

Fig. 4.— Spectrum of WISE 1118+31 obtained 
with TripleSpec (black) and LUCIFER (red). Re- 
gions of strong telluric absorption are shown in 
grey. Prominent absorption bands of CH4 and 
H 2 are indicated. The agreement between the 
two spectra is excellent. 



( Mason et al.l Il995l) known prior to the discov- 
ery of WISE 1118+31. Visible to the unaided 
eye a short distance from the Big Dipper, this 
telescopic double star was among the first to be 
recog nized as a grav itationally bound binary sys- 
tem (Hcrschcl J1804 ) . Given the 60 year period of 
the vis u al pai r it was not until 23 years later that 
Struvd (J18271 ) calculated a formal orbit. Subse- 
quently bot h components of th e visual pair with 



[Fe/H] 



-0.32 ± 0.05 dex. 



spectroscopic binaries. IHeinta ([19671 ) finds peri- 



995h 

inta 



2.536" (Maso_n_et_aLl 1995) were found to be 



ods of 669.1 days for the Aa sy s tem a nd 3.9805 
days for the Bb system. iGriffinl ( 19981 ) gives ve- 
locity amplitudes, Ka<i — 4.85 ± 0.14 km/sec and 
Ksb = 4.33 ± 0.09 km/sec in the 60 year orbit 
which lead to a dynamical parallax estimate of 
0.126" ± 0.0 023". 

However ISoderhiehrj (|l999l ) re-analyzed the 
Hipparcos data combined with other data for vi- 
sual binaries and gives 0.1197" ± 0.0008" for the 
parallax of £ UMa, and total mass for the AaBb 
system of 2.62 M . 

The weighted mean of the dynamical parallax 
and the Hipparcos parallax is 0.1206" ± 0.00074" 
but with x 2 — 8.6 for 1 degree of freedom in the fit 
for the mean, so we inflate the errors by a factor of 
V8£ and adopt 0.1206"±0.0022" for the parallax. 

Since the poss ible 5 th component, seen by 



Mason et al.l (J1995I ) separated by 56 milli-arcseconds 
from the B component of the visual binary, was 
only detected at one epoch we will not consider it 
to be a member of the system. 

iBakos et al.l (|2002j ) give a proper motion of 
0.75"/yr in position angle 217.49°, based on 
a 29 year interval. These values resolve into 
-0.456"/yr in RA and -0.595"/yr in Dec which 
we adopt. 

3.1.1. Spectral Types and Metallicities of the 
Central Stars 



Cavrel de Strobel et al.l (|1994j ) used high reso- 
lution spectroscopy to find effective temperatures 
and gravities of 5950 ± 30 K, \ogg = 4.3 ± 0.2 
for £ UMa A; and 5650 ± 50 K, logg = 4.5 ± 0.2 
for B. These temperatures and gravities are con- 
sistent with spectral types of F 8.5V and G2V as- 
signed to £ UMa A and B by Keenan fc McNeill 
(|l989l ). ICavrel de Strobel et al.l (|l994h found that 
both stars had slightly sub-solar iron abundances, 



3.1.2. Chromo spheric Activity and Age of the 
Central Stars 



Cavrel de Strobel et al.l (|l994l) see chromo 



spheric calcium emission lines for £ UMa B but 
not for £ UMa A. They propose that the short 4 
day period orbit of the Bb system is driving the 
chromosp heric activity in B. From the lac k of emis- 
sion in A. ICavrel de Strobel et al.l (|l994l ) conclude 
that £ UMa is older than stars i n the cluste r NGC 
752, which has an age of 2 Gyr. JBall et al.l (|2005h 
report X-ray observations of £ UMa which show 
that all of the observed X-ray emission is coming 
from the 4 day period binary Bb. X-ray flux from 
component A is more than 2 orders of magnitude 
fainter, implying logLx < 27.5. Since component 
A is close to solar luminosity the ratio of X-ray 
to bolometric flux is Rx = \og(Lx/Lt, i) < —6 
which gives an age greater than 4 Gy r usin g 
Equation(A3) in iMamaiek fc Hillenbrandl (|2008l ). 
The low space velocity (U, V, W) = (-2.36 ± 
0.33, -28.45 ± 1.22,-20.26 ± 0.2 5) km/sec cal- 
culated by iKaratas et al.l (J2004T ) implies that 
£ UMa is a member of the thin disk, so ages 
much greater than 8 Gyr ar e unlikely. Finally 
ICavrel de Strobel et al.l |l994l ) found that the lu- 
minosities and effective temperatures were con- 
sistent within the errors with the masses derived 



by IHeinta (|l967l ) for a 5 Gyr isochrone calculated 
with sub-solar abundance. 

3.2. Astrometry of the Companion 

Astrometric information was extracted from 
the observed images at the various epochs us- 
ing the standard maximum likelihood technique 
in which a point spread function (PSF) was fit 
to each source profile. The positional uncertain- 
ties were estimated using an error model which 
includes the effects of instrumental and sky back- 
ground noise and PSF uncertainty. The PSF and 
its associated uncertainty map were estimated for 
each image individually using a set of bright stars 
in the field. In order to minimize systematic ef- 
fects, our astrometry was based on relative posi- 
tions, using as a reference the nearby star 2MASS 
J11183876+3125441 at a separation of approxi- 
mately 6.5" (star 1 in Figure[2]). This is sufficiently 
close that systematic errors due to such effects as 
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Fig. 5. — Astrometric data and fits for WISE 
1118+31. The blue curves and point are for 
ground-based or Low Earth Orbit observatories, 
while the red curves and points are for Spitzer. 
The bold dashed lines show the fit with proper 
motion and parallax as free parameters, while the 
lighter solid curves show the fit forced to match £ 
UMa. These fits are very similar. Note that AS 
has been displaced by a constant for clarity. 

focal-plane distortion, plate scale, and rotation er- 
rors cancel out in the relative position to an ac- 
curacy much greater than the random estimation 
errors. We do, however, assume that the proper 
motion and parallax of the reference star are negli- 
gible compared to those of WISE 1118+31. In sup- 
port of this assumption, we find no significant dif- 
ference in relative motion of WISE 1118+31 when 
the differential astrometry is repeated using more 
distant 2MASS reference stars, at separations of 
69" and 105". The astrometric esti mation proce- 
dure i s discussed in more detail by iMarsh et al 
(|2012h . 



Table [5] gives the measured position offsets 
which were input to a parallax and proper motion 
code that handles both Earth-based and Spitzer 
observations. Table |3] gives the output from the 
code for three different sets of free parameters. 
In one case the proper motion and parallax were 
forced to zero. In the second case they were forced 
to match £ UMa. In the final case the proper mo- 
tion and parallax were left as free parameters. 

The astrometric fits for WISE 1118+31 are a 
very good match to the motion of £ UMa. The 
Ax 2 = 241.0 between the best fit and a fit with 
a fixed position shows that the motion of WISE 
1118+31 has been detected with a SNR of 15.5, 



while the A% 2 = 1.28 for 3 extra degrees of 
freedom between the best fit and a fit forced to 
match the proper motion and parallax of £ UMa 
is perfectly consistent with WISE 1118+31 being 
a bound member of the £ UMa system. Figure [5] 
shows two fits: one forced to match £ UMa and 
one with the proper motion and parallax as free 
parameters. 

3.3. The Nature of the Companion 

3.3.1. Spectral Classification 

As shown in Figure 21 the spectrum of WISE 
1118+31 exhibits deep absorption bands of CU4 
and H2O indicative of T dwarfs. We derive a 
more precise spectral type using the spectral clas- 



sification schemes of iBureasser et all (12006!) and 



;pc< 



the ex tension to this system by Cushing et al 



(120111 ) whereby UGP S J072227.51-0540 31.2 (here- 
after UGPS 0722-05. lLucas et alJlJblbT) is defined 
as the T9 spectral standard and WISE 1738+2732 
is defined as the Y0 spectral standard. As shown 
in Figure El WISE 1118+31 has a spectral type 
of T8.5 based on the width of the J-band peak at 
1.27 jum. 

3.3.2. Absolute magnitude, Luminosity and Mass 

For our adopted parallax the absolute magni- 
tude is M W2 = W2 + 51og(10ro) = 13.715 ± 
0.051. This agr ees very well with the M wi vs. 
type relation in iKirkpatrick et al.l ( 20111 ) which 
gives M\Y2 = 13.71 ± 0.21 for a spectral type 
o f T8.5 ± 0.5. T he M W2 vs. type relation in 



Kirkpatrick et al.l (120121) with W1828+2650 ex- 
cluded gives M\Y2 = 13.64+0.31 which also agrees 
with the observed W2 flux. This M W2 = 13.715 
implies vL v — 1.14 x 10~ 6 L G at 4.6 /im. 

To convert vL v to L we need the factor 
vF v {W2) I '-Fbob which is equivalent to the bolo- 
metric correction, s ince vF v — (Ai SO F£// c )10 _a4 ™ 1 
((Wright et al.l2010l ) which is 1.118xl0~ 7 erg/cm 2 /sec 
for W2=0 with a v~ l spectrum, and Fboi = 
2.48 x io- 5 10-°' 4mb °i erg/cm 2 /sec, so 

2.5 log f vF ^ W2 ^ ) =-5.87 + TO b oi-m(IU2). 

V ^bol / 



(1) 



Thus 



BC{W2) 



rriboi - m(W2) 



Table 2: Astrometric Observations of WISE 1118+31. 



Date 


Aa cos S ["] 


AS ["] 


Observatory 


Band 


2010.39 


-0.306 ±0.162 


-0.444 ±0.182 


WISE 


W 


2010.91 


-0.165 ±0.211 


-0.603 ±0.191 


FanMt 


Y 


2010.91 


-0.208 ±0.162 


-0.477 ±0.184 


WISE 


W 


2011.05 


-0.232 ±0.067 


-0.794 ±0.067 


Spitzer 


ch2 


2011.12 


-0.414 ±0.234 


-0.600 ±0.194 


FanMt 


J 


2011.24 


-0.393 ±0.113 


-0.765 ±0.100 


FanMt 


J 


2011.38 


-0.651 ±0.086 


-0.829 ±0.084 


MtBglw 


J 


2011.38 


-0.554 ±0.194 


-0.723 ±0.190 


MtBglw 


H 


2011.85 


-0.639 ±0.190 


-1.280 ±0.170 


FanMt 


J 


2012.00 


-0.650 ±0.005 


-1.253 ±0.004 


WIYN 


J 


2012.00 


-0.676 ±0.010 


-1.235 ±0.012 


WIYN 


H 


2012.00 


-0.766 ±0.037 


-1.239 ±0.038 


WIYN 


K 


2012.20 


-0.850 ±0.165 


-1.339± 0.153 


FanMt 


Y 


2012.15 


-0.695 ±0.117 


-1.631 ±0.126 


Spitzer 


ch2 


2012.52 


-1.081 ±0.136 


-1.611 ±0.139 


Spitzer 


ch2 



Note.— Offsets relative to ll' 1 18 m 38.69 s +31° 25' 37.7" (J2000). 
Reference star: U h 18 m 38.77 s +31° 25' 44.2" (2MASS). 



Table 3: Astrometric Fits for WISE 1118+31. 



Type 



X 



#df cos Sda/dt ["/yr] dS/dt ["/yr] 



Fixed 270.836 28 

Forced to £ UMa 29.856 28 -0.456 -0.595 0.1206 

Free 28.579 25 -0.419 ± 0.048 -0.563 ±0.045 0.124 ±0.033 



2.5 log 



/ vF„(W2) \ 



V F, 



bol 






5.87. (2) 



Note that vF v /F\> \ is non- negative and also 
normalized since §{yF v /F\ 10 \)d\a.i> — 1 by defi- 
nition. Since the normalization integral is domi- 
nated by the brightest peaks, errors in the faint 
parts of the spectral energy distribution (SED) 
have only small effects on the derived value of 
vF u (W2)/Fb i or equivalently the bolometric cor- 
rection. The brightest peaks in vF v jF\, Q \ for late 
T dwarfs are in the J band and the W2 band, so 
if we get the J — W2 color right then our value of 
vF l/ (W2)/Fb i will be fairly accurate. 

We will use a model spectrum to evaluate 
vF v {W2)/F\> \ in order to fill in the gaps between 
the photometric passbands, but we are only us- 
ing the shape of the spectrum to calculate the 



bolometric correction. Thus, even if the effective 
temperature T e g is a parameter of the model, that 
does not imply that the model parameter is the 
actual T e g of the star. 

Given a model spectrum F v , the ratio vF v (W2)/ 'Fboi 
is given by 



vF u 



${vF v /hv)R(v)dh\v 



F 



boi J ([1 erg/cm ]/hv)R{v)dh\v J F v dv 
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wher e R{y) is the response per photon (jWright et al. 
20101 ). The model we use is a linear combination 



of sulfide cloud models from iMorlev et al] (J2012J) 

with T eff = 600 K, log(.g) = 5.0. We form the sum 
of 0.2 times the flux from the cloudless model and 
0.8 times the flux from the f se d — 4 model. Phys- 
ically this is a brown dwarf that is 20% covered 
by clear zones and 80% covered by cloudy bands. 



c 

CC 

H — ' 

CO 

c 
o 
O 



CM 

7 

CD 
CM 




1.3 1.4 1.5 1.6 
Wavelength (|xm) 



1.8 



Fig. 6. — Sequence of spectral standards from 
T7 to YO (red) along with the spectrum WISE 
1118+31 (black). The spectral standards are 
2MA SS 0727+1710 (T7, iBurgasser et~all fc oOfx 



2002) and 2 MASS 0415-0935 (T8, Burgasscr et al. 



2006, 



2011 



(Y0, 



2002) , UGPS 072 2-05 (T9, ICushing et al 



Lucas et al 



Cushing et al 



20101) . and WISE 1738+2732 



120111) . The spectrum of WISE 
1118+3125 has been smoothed to a resolving 
power of \/d\ = 1000 for display purposes. Spec- 
tra are normalized to unity over the 1.26-1.27 
/xm wavelength range and offset for clarity (dot- 
ted lines). Regions of high telluric absorption are 
shown in grey. WISE 1118+31 is classified as 
T8.5+0.5. The mismatch between the spectrum 
of WISE 1118+31 and the templates at -1.63 fjm 
is due to poor subtraction of the OH sky lines. 



This linear combination gives vF v jF\, \ = 1.46, 
Jmko -W2 = 4.60, (Y - J)mko = 1-23, 
(J - H) MKO = -0.18, (J - K s ) MKO = -0.06 
and Wl - W2 = 3.17. The model matches the 
observed J — W2 color by design but is too bright 
in K s and too faint in Wl. If we make an ad hoc 
correction to the SED by making the wavelength 
range from 3.0 to 3.8 fim covering the Wl band 
0.4 mag brighter that would match the Wl — W2 
color and change the normalization integral by 
2%. If we make the wavelength range from 2.0 
to 2.3 /im covering the K s band 0.8 mag fainter, 
then the J — K s color matches the data and the 
normalization integral changes by -1.3%. Thus we 
adopt vF v {W2)/F ho \ = 1.46 with a 10% uncer- 
tainty to allow for such errors in the model, giving 
L = (1.14/1.46) x 10- 6 = 10- 6107 ± 0043 L G and 
BC(W2) = 6.28 so M bol = 20.0 ± 0.11 for WISE 
11 18+31. Interpolating the L vs. M and t figures 
in lSaumon fe Marlevl ( 20081 ) gives this luminosity 
for masses of 14 to 31 Mj with ages of 2 and 8 
Gyr using the cloudy models, and 17 to 38 Mj for 
the cloudless models. 

Even though the model parameters do not nec- 
essarily apply to the actual star, we find that 
they do give a mass and luminosity consistent 
with those derived a bove. The T n ff = 6 00 K , 
log(g) = 5.0 models of iHubenv fc Burrows! (|2007h 
have radii of 0.91 Rj, and a mass M — gR 2 /G = 
32 Mj. The luminosity calculated directly from 
T cS and R is L = ^R 2 T^ S = 10" 601 L© so 
this model is reasonably self-consistent. Calcu- 
lating the temperature from R and L gives T c s — 
(0.91 J R J /i?) 1 / 2 (567±14) K. 

3.4. Another Binary? 

If something like the iMorlev et al.l ( 20121 1 sul- 



fide clouds is not causing the red J — W2 color of 
WISE 1118+31, then there is a tension between 
the absolute magnitude My/2, which requires a 
fairly high T c g, and the colors and spectral type 
which suggest a lower temperature. This could 
be ameliorated if WISE 1118+31 were an equal 
mass binary which would increase the total ra- 
diating area. In this case the absolute magni- 
tude of one component of the binary would be 
Mw2 = 14.47 and th e colors would remain the 
same. The best fitting IHubenv fc Burrows! (|2007l ) 
model is then a non-chemical equilibrium model 
with T off = 500 K and log( 5 ) = 5.0. This model 



has M W2 = 14.27, W1-W2 = 4.43, J(MK0)-W2 
= 4.57, and R = 0.888 Rj. vF v jF ho \ = 1.546 in 
the W2 band for this model, so the luminosity of a 
single component is L = 0.5 x 1.14 x 10~ 6 /1.546 = 
3.7 x 10 -7 Lq. The mass range fo r ages from 2 
to 8 Gyr in the ISaumon fe Marlevl (J2008I ) cloud- 
less models is 12 to 29 Mj, while the cloudy mod- 
els give 10 to 25 Mj. Thus the range of total 
masses is 21 to 59 Mj under the equal mass binary 
hypothesis, which is not much different than the 
range under the single object hypothesis. How- 
ever, the absolute magnitude of the components 
would no longer agree with the Mw 2 vs. type re- 
lation (JKirkpatrick et al.ll20lil l2012h . 

As stated in JJ2.2.51 we saw no evidence for bi- 
narity in the 14 Apr 2010 NIRC2 images of WISE 
1118+31. To determine the upper limit on the 
separation of an equal brightness companion, we 
added a shifted copy of the mosaic to the original 
mosaic. At separations < 4 pixels, we are unable 
to convincingly establish the presence of the com- 
panion because of the difficulty in differentiating 
between an irregular PSF caused by an off-axis tip- 
tilt reference star and an elongated PSF caused by 
a companion when no other stars are present on 
the array. At separations of >5 pixels (>50 mas), 
the presence of a companion is obvious. Therefore, 
we set the upper limit on the separation to 50 mas, 
slightly less than the FWHM of the image. If this 
was not due to a temporary conjunction, so the 
actual separation were less than 0.4 AU, then the 
orbital period could be less than 2 years. 

3.5. The View from the Neighborhood of 
WISE 1118+31 

With a projected separation of 4100 AU from 
the four main components of the £ UMa system, 
WISE 1118+31 has a distant but interesting per- 
spective. The two binaries that form the 2" visual 
pair as seen from the Earth are separated by 20 
AU, or by an angular scale of 15' from the distant 
perspective of WISE 1118+31. At that distance 
the A and B components would each shine with 
an apparent visual magnitude of -9, one hundred 
times brighter than Venus in the skies of Earth. 

Given the complexity of the £ UMa system, it is 
reasonable to speculate that WISE 1118+31 could 
be a component that was once more closely bound 
to the system, which was thrown into an orbit with 
a very large apocenter during a three body inter- 



action that tightened the £ UMa Bb binary. This 
would make WISE 1118+31 similar in history to 
an Oort Cloud comet, formed at a radius of tens 
of AU from the central star system, lifted into an 
orbit with a high apocenter by multi-body per- 
turbations, fol]owed_bygalactic tides raising the 
pericenter ( Duncan et al.lll987l ). Given an orbital 
period of order 10 5 years, the orbital eccentricity 
of WISE 1118+31 may ultimately be detectable. 
With an apparent orbital radius of 8.5' the orbital 
motion will amount to tens of milliarcseconds per 
year relative to the primary system, and a radial 
velocity difference of about 1 km/sec. Characteriz- 
ing this motion, specifically the differential proper 
motion between WISE 1118+31 and its primary, 
is tractable given modern infrared and visual as- 
trometric capability, and certainly will be accom- 
plished with the passage of time. 

4. Conclusions 

Alula Australis (£ UMa), a solar neighbor- 
hood visual binary where each component is it- 
self a spectroscopic binary, possesses an ultra-cool 
brown dwarf (T8.5) companion at a projected sep- 
aration of 4100 AU. This system is similar to, but 
somew hat more complex than, the binaries stud- 
ied by lAllen et all (J2012I) . who found that 1 in 
5 spectroscopic binary systems ha d distant com- 
mon proper motion companions. iFahertv et al 



([20101 ) found that wide companions were much 



more likely with binary or tertiary central objects, 
and £ UMa is a continuation of that trend. Thus 
the £ UMa system is uncommon primarily in the 
apparent magnitude of the central star system, 
due to its proximity to the Sun. The £ UMa sys- 
tem provides a very accurate absolute magnitude 
for a T8.5 brown dwarf, but both the age of the 
system and the mass of the brown dwarf are un- 
certain due to the difficulty of determining the age 
of a star in the middle of its main sequence life- 
time. WISE 1118+31 appears to be slightly less 
luminous and redder than the T8 dwarf WISEP 
J1423+01 which is a co mmon proper motion com- 



panion to BD+01 2920 (jPinfield et al.l l2012) 



The £ UMa system is amenable to detailed 
study since it is quite close to the Sun. An 
accurate (sub km/sec) radial velocity for WISE 
1118+31 would be valuable information for deter- 
mining its orbit around £ UMa. Improved proper 
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motions for both WISE 1118+31 and the center 
of mass of the central quadruple star system are 
also needed for orbit characterization. 

There is a tension between the red J-W2 color of 
WISE 1118+31 and its absolute magnitude My/% 
when fitting to older models. This discrepancy 
could be relaxed if WISE 1118+31 were a binary, 
but one AO observation showed no evidence for 
bin arity. Thus the new types of clouds considered 
by iMorlev et al.l ( 2012 ) may well be significant, 
and we conclude that model spectra need to be 
updated to fit the observations of WISE 1118+31 
and the many other late T and Y dwarfs found in 
the WISE data. 
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